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THE  ENERGY  CRITERIA  FOR  STABILITY  OF  STRUCTURES 


by  Gerald  Wempner 
Member  ASCE 


ABSTRACT 

The  energy  criteria  of  Trefftz  and  Koiter,  for  the  critical  load  upon  an  elastic  system  and  for 
the  stability  of  the  system  at  the  critical  load,  are  presented.  The  presentation  employs  geometrical 
interpretations  and  simple  examples  to  exhibit  the  essential  features  of  the  criteria  and  the  related 
behavior  of  the  structural  system. 

INTRODUCTION 


If  the  load  upon  a  structure  attains  a  critical  value,  the  structure  buckles.  The  buckling  may 
entail  a  gradual,  albeit  excessive,  deflection  as  the  load  exceeds  the  critical  value.  Otherwise,  the 
buckling  may  mean  an  abrupt  collapse,  so-called  snap-buckling.  The  former  occurs  if  an  elastic 
system  is  stable  at  the  critical  load  and  the  latter  if  it  is  unstable.  From  a  practical  viewpoint,  the 
snap-buckling  is  the  more  dangerous  phenomenon.  Moreover,  the  structure  which  exhibits 
snap-buckling  is  also  sensitive  to  imperfections;  that  is,  the  snap-buckling  of  a  real  (imperfect) 
structure  may  occur  at  loads  much  less  than  the  critical  load  of  the  ideal  structure. 

The  phenomenon  of  snap-buckling  is  especial!,  prevalent  in  thin  shells  and,  curiously,  the  most 
efficient  shell-like  structures  are  the  most  susceptible  to  such  catastrophic  failure.  Consequently,  as 
our  attention  turns  increasingly  to  thin  shells,  to  reduce  costs  and  weight  and  to  achieve  structural 
and  esthetic  aims,  the  questions  of  stability  and  imperfection  sensitivity  are  paramount. 

The  determination  of  the  critical  load  by  a  stationary  criterion  upon  the  potential  energy  was 
given  by  Trefftz*  in  1933.  However,  the  question  of  stability  at  the  critical  state  and  the  subsequent 
behavior  remained  unanswered. 

A  most  significant  work  on  the  questions  of  stability  at  the  critical  load,  post-buckliri  j 
behavior  and  imperfection  sensitivity,  was  the  thesis  of  W.  T.  Koiter^  in  194b.  By  employing  a 
variational  approach  and  the  criteria  for  a  minimum  of  the  potential  energy,  Koiter  developed 
stationary  conditions  for  stability  at  the  critical  load.  In  addition,  he  examined  the  effects  of  small 
geometrical  imperfections  and  showed  how  such  imperfection  can  drastically  reduce  the  buckling 
loads  upon  real  structures. 


Professor,  The  University  of  Alabama  in  Huntsville 


>  “’i  '  V'SC-*-, 


tsst* 


Although  Koiter’s  thesis  is  now  available  in  an  English  translation^,  the  rigorous  character  of 
the  work  seems  to  inhibit  a  widespread  appreciation  and  usage.  By  appealing  to  various  geometrical 
interpretatic  is,  the  following  presentation  offers  a  simple  introduction  to  the  essential  features  and 
the  consequences  of  Koiter’s  work.  For  simplicity,  the  ideas  are  developed  here  for  the  discrete 
mechanical  system,  but  are  readily  extended  to  the  continuous  system. 


PRIMCIPI.ES  OF  STATIONARY  AND  MINIMUM  POTENTIAL  ENERGY  i 

1 

If  a  discrete  conservative  mechanical  system  has  N  degrees-of-freedom,  the  configuration  is  'j 

determined  by  N  generalized  coordinates  qj  and  the  potential  energy  is  a  function  of  the 

< 

coordinates: 

>  • 

V  =  V(qj,  -,  qN).  1 

Throughout  our  development,  the  forces  are  assumed  continuous  with  continuous  derivatives,  and 
so,  the  potential  is  also  assumed  continuous  with  continuous  derivatives. 

The  principle  of  virtual  work  asserts  that  a  state  (qj)  is  a  state  of  equilibrium  if  the  potential  V 
(qj)  is  stationary:  that  is,  a  small  displacement  (q,  —  q,)«  produces  no  change  of  first-degree  in  the 
potential: 

W 

AVi%T(qi"qi)  =  0  2 

oq, 

where  the  bar  (— )  signifies  evaluation  in  the  reference  state  (qj) 


The  motionless  system  is  in  a  state  of  stable  equilibrium  if  and  only  if  the  potential  V  is  a 
proper  minimum^;  that  is,  the  state  (qj)  is  a  state  of  stable  equilibrium  if,  and  only  if, 


AV  =  V(q1(--,  qjM)  ~  V(qj,  -,  qN)  >  0  3 

for  all  displacements  (q,  -  q,).  sufficiently  small.  The  qualification,  sufficiently  small,  is  added  to 
emphasize  that  we  require  only  the  local  minimum.  For  example,  a  ball  resting  in  a  shallow  valley  is 
in  a  stable  configuration,  strictly  speaking,  but  a  small  jolt  may  kick  the  ball  over  the  adjoining  hill 
and  into  a  lower  valley  (a  more  stable  position). 

AN  EXAMPLE  OF  STRUCTURAL  INSTABILITY 


The  system  depicted  in  Fig.  1  is  composed  of  two  rigid  links  AB  and  BC,  joined  by  a 
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frictionless  pin  at  B  and  constrained  by  a  linear  extensional  spring  k  and  torsional  spring  0.  The 
extensianal  spring  resists  lateral  displacement  W  with  a  force  F  =  kW  and  the  torsional  spring  resists 
the  relative  rotation  (20)  with  a  couple  C  =  0(2 0).  The  top  A  is  constrained  to  move  vertically  while 
the  bottom  C  is  pinned  to  a  fi:;ed  support.  Consider  the  equilibrium  of  this  system  under  the  action 
of  an  axial  force,  P  =  constant,  applied  to  the  end  A: 


Our  system  has  one-degree-of-freedom.  The  configuration  is  determined  by  the  kinematic 
variable  0  and  the  total  potential  energy  of  the  system  is 

kL^  o  9 

V  =  -  sin'fy  +  200*-  +  2PLcos0  4 

2 

By  the  principle  of  stationary  potential  energy,  the  system  is  in  equilibrium  if 

klJ*  sin  0  cos 0  +  4 00  —  2PL  sin0  =0  5 
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Evidently,  0=0  determines  an  equilibrium  configuration  as  it  satisfies  (5)  for  all  choicest,  L, 
0,  P.  Now,  we  ask:  Are  there  other  solutions  Of  0  which  satisfy  (5)?  If  so,  we  can  divide  (5)  by  2L 
sin0  and  obtain  , 


P  = - cos0  + -  "  6 

2  L  sin0 

Some  plots  of  (6)  trace  the  solid  lines  in  Fig.  2.  The  ordinate  is  the  dimensionless  load  P/Pcr  where 


kL  2 0 
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At  the  point  labeled  0  in  Fig.  2  there  is  a  bifurcation;  one  branch  is  the  vertical  line  0=0  and 
another  branch,  is  the  curve  of  (6).  i  : 

Let  us  consider  th^  system  when  kL*  f  8/3  and  examine  its  behavior-as  the  load  P  is  gradually 
'applied.  The  'load-deflection  curve  is  the  solid  curve  'of  Fig,  3;  the  dotted  curves  are 
energy -deflection  curves  for  P  =  Pcr  and  P  =  0.8  Pcr  : 


oi'  \ 


POTENTIAL 


[OJ 


The  system  can  and  does  sustain  the  load  in  the  straight  configuration  0=0  until  the  load  reaches 
the  critical  value  P  =  Pcr  For  loads  P  ^  Pcr  the  configurations  0  =  0  are  unstable.  At  the  critical 
state  0,  a  slight  additional  load  will  cause  the  system  to  snap-thru  to  the  configuration  of  point  Q, 
because  the  total  potential  energy  at  Q  is  less  than  at  0;  point  Q'lies  below  point  O'.  Moreover, 
every  intermediate  configuration  has  a  lower  energy  level.  In  particular,  the  configurations  adjacent 
to  0  =  0  have  less  potential  energy  and,  consequently,  the  system  tends  to  move  from  the  straight 
configuration.  This  is  apparent  from  the  energy  curve  near  O'in  Fig.  3. 

Now  consider  the  same  system  as  the  load  is  gradually  applied  until  P  =  0.8  Pcr  The  system 
sustains  this  load  in  the  straight  configuration  of  point  A  in  Fig.  3.  At  this  load  the  potential  energy 
follows  the  curve  A'B'C'  with  relative  minima  at  A'  and  C'  and  a  relative  maximum  at  B'.  The 
configurations  of  A  and  C  are  stable  while  that  of  B  is  unstable.  If  some  energy  were  supplied  to  the 
system  in  configuration  A,  the  system  could  be  kicked  over  the  energy  hill  A'B’  into  the  valley  at  C' 
corresponding  to  C  on  the  load-curve.  The  configuration  of  A  is  less  stable  than  that  of  C  because  a 
slight  disturbance  can  cause  a  violent  snap-through  from  A  to  C. 

Notice  that  the  curve  OS  in  Fig.  2  has  a  positive  slope  everywhere.  It  represents  stable 
configurations.  However,  because  the  deflections  increase  rapidly  with  load,  a  structure  of  this  kind 
may  be  unusable  for  loads  beyond  the  critical  value,  i.e.  P  •*  Pcr  A  column,  for  example,  is  said  to 
buckle  when  the  load  exceeds  the  critical  value  of  the  bifurcation  point. 

At  the  critical  load  the  unbuckled  configuration  ( 0  =  0)  may  be  stable  or  unstable;  that  is,  the 
system  may  sustain  additional  load  accompanied  by  a  gradual  increase  in  deflection  (curve  OS  of 
Fig.  2)  or  it  may  snap  abruptly  to  a  severely  deformed  configuration  at  the  slightest  disturbance 
(curve  OQ  of  Fig.  2).  From  a  practical  viewpoint  the  question  of  stability  at  the  critical  load  is 
extremely  important.  In  the  present  axamplc,  the  question  is  easily  resolved  by  examining  the 
load-deflection  function  P  =  P{0),  However,  we  are  expressly  concerned  with  the  energy  criteria  and 
so  we  examine  the  conditions  for  a  minimum  of  potential  energy. 

In  the  neighborhood  of  the  reference  state,  the  potential  V  of  (4)  can  be  represented  oy  the 
series  expansion: 

kL2  o  1  PL  o  a 

V  =  ( - +  2/1  ~  PL)0  z+  —  ( - kL2)04  +  --  bu 

2  6  2 

=  v2  +  v4  +■*  8b 


6 


«T  tvp^iuwgfmyamrvmi  »»wrg»?WftSE»?!S!5SSSSS>  *  nf  -yo— 


where  V^  denotes  the  term  of  degree  N  in  the  variable  0 .  Notice  that  the  odd  powers  are  absent 
because  the  structure  is  symmetrical;  equal  deflections  to  the  right  or  left  produce  the  same  change 
of  potential.  If  V 2  0,  then  V2  dominates  and,  sufficiently  near  to  the  reference  state  0=0,  the 

stationary  condition  (5)  can  be  replaced  by 

d(Vo)  kL2 

_L_2i  =  2{ - +  2/3  —  PL)  0=0  9 

dO  2 

The  ‘‘equilibrium"  condition  (9)  has  a  nontrivial  solution  Of  0,  if  and  only  if  the  parenthetical 
factor  vanishes,  that  is,  if  the  load  has  the  critical  value  (7).  The  stationary  condition  upon  the 
second-decree  term  Vp,  is  the  criterion  of  Trefft2. 


At  the  critical  state,  V 2  “  0  and 


V  -  v4  + 


The  system  is  stable  according  to  the  principle  of  minimum  energy,  if 


v4  -  0 


In  accordance  with  (7)  and  (8)  the  system  is  stable  if 


If  the  system  of  Fig.  1  is  imperfect,  say  the  linkage  has  an  initial  angle  0Q,  then  a  lateral 
deflection  W  accompanies  the  initial  increment  of  load.  As  the  load  increases,  a  plot  of  load  versus 
deflection  traces  the  dotted  curve  of  Fig.  2  and  approaches  the  solid  trace  of  the  perfect  system.  If 
the  system  exhibits  the  snap-through  characteristic,  then  the  load-deflection  curve  has  a  negative 
slope,  as  OB,  and  the  crest  of  the  actual  (dotted)  curve  falls  below  the  line  P  =  Pcr.  This  suggests 
that  structures  exhibiting  the  snao-throuoh  phenomena  are  also  sensitive  to  imperfections. 
Experimental  evidence  confirms  our  suspicions. 

STABILITY  OF  A  DISCRETE  MECHANICAL  SYSTEM 


The  essential  features  of  the  criterion1  for  a  critical  load  and  the  criteria'1  for  stability  at  the 
critical  load  are  exhibited  most  clearly  by  a  discrete  mechanical  system.  The  underlying  concepts 
apply  to  a  continuous  system  so  that  the  criteria  are  readily  extended.  The  characteristics  and 
behavior  of  our  system  follow: 


7 


T/WfTT^r^i*'  f  vTJFrV3s  Wi  OT  IV/?  ^  v  tv  &gV£j  V 

***"?*» « **pg»«$7gy r  T^*iy*<yy^ 


All  loads  upon  the  system  are  assumed  to  increase  in  proportion  and,  therefore,  the  magnitude 
is  given  by  a  positive  parameter  X.  A  configuration  of  the  system  is  defined  by  N  generalized 
coordinates  q j (i  =  1,  --,N).  As  the  loading  parameter  is  increased  from  zero,  the  equilibrium  states 
trace  a  path  in  a  configuration-load  space.  For  example,  a  system  with  two  degrees  of  freedom  (q  j, 
q2>  follows  a  path  in  the  space  (q^,  q2^)  of  Fig.  4a  or  Fig.  4b. 

The  point  P  of  Fig.  4a  or  Fig.  4b  is  a  critical  state  characterized  by  the  existence  of 
neighboring  states  which  are  not  uniquely  determined  by  an  increment  of  the  load.  At  the  critical 
point  P  of  Fig.  4a  the  path  OP  forms  two  branches,  PR  and  PQ.  The  branch  PQ  may  ascend  or 

A 

descend,  cf  the  tangent  W  may  be  normal  to  the  X  axis.  At  the  point  P  of  Fig.  4b,  the  smooth  curve 
has  a  ta^ant  W  normal  to  the  X  axis.  Then,  if  t  denotes  arc  length  along  the  path,  the  path  PQ  at  P 
is  characterized  by  the  condition  dX/dc  =  0;  in  words,  the  system  tends  to  move  from  P  with  no 
increase  of  load. 

The  critical  state  of  Fig.  4a  occurs  at  a  bifurcation  point;  two  paths  of  equilibrium  emanate  from  P. 
However,  the  path  PR  represents  unstable  paths  which  can  not  be  realized.  Actually,  the  system 
tends  to  move  along  PQ.  If  PQ  is  an  ascending  path,  then  additional  loading  is  needed,  and  the 
system  is  said  to  be  stable  at  the  critical  state.  In  actuality,  a  very  slight  increment  is  usually  enough 
to  cause  an  unacceptable  deflection  and  the  system  is  said  to  buckle.  If  the  curve  PQ  is  descending, 
the  system  collapses  under  the  critical  load  X*. 

The  path  of  Fig.  4b  is  entirely  smooth,  but  reaches  a  so-called  limit  point  P.  The  state  of  P  is 

A 

again  critical  in  the  sense  that  the  tangent  W  is  normal  to  the  X  axis.  At  P  the  system  tends  to  move 
under  the  critical  load  X*.  It  tends  to  buckle,  but  it  is  theoretically  stable  if  PQ  is  an  ascending 
curve.  It  collapses  if  the  path  PQ  descends. 

By  our  remarks,  instability  is  signaled  by  the  advent  of  excessive  deflections  which  are 
produced  by  a  critical  load  X*.  However,  the  stability  of  a  conservative  system  can  be  characterized 
by  an  energy  criterion:  The  conservative  mechanical  system  is  in  stable  equilibrium  if  the  potential 
energy  is  a  proper  min, mum,  unstable  if  any  adjacent  state  has  a  lower  potential.  Let  us  apply  the 
energy  criterion  at  the  critical  state: 


TREFFTZ  CRITERION  FOR  THE  CRITICAL  STATE 


We  presume  that  the  potential  energy  can  be  expanded  in  a  power  series  about  the  critical 
state.  If  (qjiX)  defines  a  state  of  equilibrium,  Uj  =  Aqj  defines  a  displacement  from  the  reference 
state,  and  V  the  change  of  potential  caused  by  the  displacement,  then 


V  -  AjUj  <•  2"Aijuiuj+3;  Aijkuiuju«  !  ”  A'jKl  uiuJuKul  *  - 
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Since  the  state  is  a  state  of  equilibrium,  in  accordance  with  the  principle  of  virtual  work, 


Aj  -  0 
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If  the  quadratic  term  of  (11)  does  not  vanish  identically,  then  it  dominates  for  small  enough 
displacement.  It  follows  that  the  state  is  stable  if 


v2(Uj)=  “AjjUjUj  >0 
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The  state  is  critical  if 


-AijUjUj.-Q 
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In  words,  the  state  is  critical,  if  there  exists  one  (oi  more),  non-zero  displacement(s)  Uj  which  causes 
the  quadratic  term  to  vanish,  i.o. 


v2(Cj)  "  “AjjUjUj  -  0 
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The  displacement  u  is  a  buckling  mode. 

A  minimum  is  characterized  by  a  stationary  condition.  Hoie,  the  requited  minimum  of  V2(Uj) 
is  determined  by  the  stationary  criterion  of  Trofftz:  For  an  arbitrary  vniiation  <Siij, 

f>V2  '  Ajjiijfuij  "0  17 

It  follows  that  the  buckling  mode  Uj  is  a  nontrivial  solution  of  the  equations: 

AjjUj  -  0  18 

The  homogeneous  system  has  a  non-tiivial  solution,  if  and  only  if  the  dotuiminant  of  coelfiuenl 
vanishes 


I  AjjtX)  1-0 

The  least  root  of  (19)  determines  the  critical  load 


19 


10 


Let 


Uj  =  eWj 


where  W:  are  components  of  the  unit  vector  in  our  N-dimensional  space  of  q;,  i.e. 


WjWj  =  1 


The  parameter  c  measures  the  magnitude  of  an  excursion  from  the  critical  state  and  Wj  defines  the 
direction  of  the  buckling. 

DIFFERENTIAL  GEOMETRY  OF  A  PATH  FROM  THE  CRITICAL  STATE 

Let  us  consider  a  movement  along  the  path  emanating  from  the  critical  point  P.  In  the  plane  of 
(qj,  q2),  we  see  a  path  as  shown  in  Fig.  5. 


In  Fig.  5,  p  denotes  the  curvature  of  the  path  PQ  at  P,  W  is  the  unit  tangent  at  P  and  V  the  unit 

/\  A 

normal.  The  displacement  from  P  to  Q  can  be  expressed  in  the  form  u  =  £W  +  j?V  or,  if  e  denotes 
arc-length  along  PQ, 

du  1  2d2u 
u  =  c — +  — e  - +  ---•- 


oc  2  dc2 


/v  C-  'y 

=  c-W  +  V  + 

2  p 


If  we  accept  an  approximation  of  second-degree  in  the  arc-length  c  then 


«  -  ;  -  -■ 


22a, b 


In  the  N-dimensional  space,  as  in  the  2-dimensional  space,  one  can  define  an  arc-length  t  along 
a  path  stemming  from  the  critical  state,  i.e. 


dUjdUj  =  dc‘ 


A  component  of  the  unit  tangent  is 


Wj  - 
1  dr- 


11 


A  component  of  the  unit  normal  is 


KVj  = 


d^uj 

de2 


24 


The  displacement  along  a  small  segment  is 


U:  = 


duj 
—  c.  +  - 

de  2 


1  d^uj 
de^ 


1  e2  + 


25 


_  e2 

cW:  + - K  V:  +  ... 

I  2  I 


Here  Vj  is  normalized  in  the  manner  of  (21). 


KQITER'S  CRITERIA  FOR  STABILITY  AT  THE  CRITICAL  LOAD 


A  small  displacement  from  the  critical  sta.e  is  given  by  (25).  The  buckling  mode  Wj  is 
determined  according  to  (18),  (19)  and  (20).  Now,  we  seek  the  normal  Vj  and  curvature  k  which 


2 


i^4 

\  >'"r*  -  ’■»*"•’■  ‘ 
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determine  the  curved  path  of  minimum  change  V.  The  change  of  potential  follows  from  (ll)jand 
(25)  and  simplifies  according  to  (13),  (16)  and  (18):  ■ 


V  =  iTAijKWiWjWK  +  7K  AijViVi  +  7AijKlWiWiWKWl  +  J*  AijKWiWiVK  +  0^>  26 


If  Vj  3  0  and  e  is  sufficiently  small,  the  initial  term  of  (26)  is  dominant.  Since  the  sign  of  the 
initial  (cubic)  term  can  be  positive  or  negative,  depending  on  the  sense  of  the  displacement  Wj,  a 
necessary  condition  for  stability  follows: 


A3  3,  AijKWiWiWK  =  0  27 

If  Ao  vanishes,  as  it  usually  does  in  the  case  of  a  symmetrical  structure,  then  the  sign  of  V  rests  j 
J  > 
with  the  terms  of  higher  degree.  If  V  is  negative  for  one  displacement  V,  then  the  system  is  j 

unstable.  The  minimum  of  (26)  is  stationary,  i.e.  >S(V)  =  0,  for  vaiiations  of  Vj.  The  stationary  * 

conditions  follow  } 


W-  Ai jK  WiWK  --+°W  28  j 

t 

If  the  terms  of  higher  degree  ate  neglected,  then  equations  (28)  constitute  a  linear  system  in  the/ 

_  : 

displacement  V  .  In  accordance  with  (23)  and  (24),  the  solution  V  is  to  satisfy  the  orthogonality 
>  '  > 

condition:  ; 


VjWj  --  0 


If  follows  from  (28)  that 


*2Vivj-  -«Aii,W«  ■  30 

-  e2- 

The  potential  change  cot  responding  to  the  displacement  u(  -  rW(  +  IS  °6tai'iod  horn  (26)  and 
simplified  by  means  (27)  and  (30): 


V  -  <  % 


where 


A4  » ,  AijK  I  WjWjjWK  W, 


8  V-Vi 
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/'  '  f  ■  *  V*'  '  V’SrW'Wi 


'"■?'ii^vC3Si:r#'f^ *“’•""  *"  vjrr 


'  r:, -2^2S7?327£~'  to.*^('*w 
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The  system  is  stable  if 


A4  >0 


The  system  is  unstable  if 


A4  <0 


In  a  system  with  one  degree  of  freedom,  Vj  =  0  and  the  final  term  of  (32)  vanishes. 


EQUILIBRIUM  STATES  NEAR  THE  CRITICAL  LOAD 


In  our  preceding  view  of  stability  at  the  critical  load  X*.  we  examined  the  energy  increment 
upon  excursions  from  the  critical  state,  but  assumed  that  the  load  remained  constant.  Such 
excursions  follow  the  path  of  minimum  potential  on  a  hyperplane  (X  =  X*)  in  the  configuration-load 
space  (qj ;X).  To  trace  a  path  of  equilibrium  from  the  critical  state  requires,  in  general,  a  change  in 
the  load.  Let  us  now  explore  states  of  equilibrium  near  the  reference  state  of  equilibrium  (qpX*). 
To  this  end,  we  assume  that  the  potential  V(qjjX)  can  be  expanded  in  a  Taylor’s  series  in  the  load  X, 
as  well  as  the  displacement  Uj.  Then,  in  place  of  (11),  we  have 

V  -  (AjUj  +  — AjjUjUj  "*■  31  AjjK ujujuK  •••)  +  (AjUj  +  — AjjUjUj  +^AjjKUjUjUK  +  ...)  (X  ~X*)  +  ...  34 

Here  the  prime  signifies  a  derivative  with  respect  to  the  parameter  X  and  each  of  the  coefficients 
(Aj,  A-,  etc.)  is  evaluated  at  the  critical  load. 

Along  a  smooth  path  from  the  reference  state  in  the  configuration-load  space,  the 
"displacement"  includes  a  component  in  the  direction  of  X,  as  well  as  the  direction  of  qj.  In  place  of 
(25),  we  have 


Uj  =  eUj  +  yicVj  +  ... 


(X  —X*)  =  eX  '  +  ~K(i  +  ... 


Here,  the  vector  (u';  X*)  is  the  unit  tangent  and  (Vj ;n)  is  the  principal  normal  at  (qj ;  X*)  of  the  path 
which  traces  equilibrium  states  in  the  space  of  configuration-load  (qj;  X). 


I  i 


Upon  substituting  (  35a, b)  into  {  34),  we  obtain 


V  =  e(AjUj)  +  e2H  AyUjUj  +  AjupO  + 


36 


The  principle  of  stationary  potential  energy  gives  the  equations  of  equilibrium  at  the  reference 
state: 


A)  =0 


37 


In  view  of  (37),  the  quadratic  terms  (c2)  dominate  (36).  The  stationary  principle  S(V)  =  0, 
gives  the  equilibrium  equations  for  states  very  near  the  reference  state: 


Ajjuj  -  -A-X' 


38 


Now,  the  reference  state  is  critical,  if 


A-X'  =  0 


39 


In  words,  either  -  0,  which  implies  the  existence  of  an  adjacent  state  at  the  same  levei  of  loading, 
and/or  A|  =  0,  which  holds,  as  (37),  if  the  reference  configuration  is  an  equilibrium  configuration 
for  X  f  X*.  Then,  the  equilibrium  equations  of  the  neighboring  state  follow: 


AjjUj  -  0 


40 


Equations  (40)  are  the  equations  (18)  of  the  Trofftz  condition  (17).  The  solution  of  (40)  is  the 
buckling  mode 


ui  =  Wj 


41 


Suppose,  firstly,  that  A^  =  0  in  (39)  and  X' /  0.  Then,  according  to  (37),  (40)  and  (41),  the 
potential  of  (34)  and  (36)  takes  the  form: 


V  =  e3(A3  +  a'2X'  +  ...)  +  0(e4) 


42 
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where 


A3  =5  AijKWiWiWK 

a2  '  J  A-jW-.Wj 


sgasrassB®H3K 


We  accept  the  indicated  terms  of  (42)  as  our  approximation  and,  therefore,  set 


eX'  =  X  -X* 


Our  approximation  of  (  42)  follows: 


V  =  e3A3  +  c2A 2<X  —  X*) 


The  principle  of  stationary  potential  provides  the  equation  of  equilibrium: 


=  3e2A3  +  2eA2(X-X*)  =  0 


(X  -X*) 


The  state  is  stable  if  the  potential  is  a  minimum,  that  is,  if 


— g-  =  6c;A3  +2A2(X~X*)  ’0 


or,  in  accordance  with  (47b),  the  system  is  stable  in  the  adjacent  state  if 


-A2  (X  -  X*)  >0 


In  accordance  with  (34),  (40)  and  (41),  the  quadratic  terms  of  V  in  the  buckled  mode  follow: 


v2(Wj)s  “  Ajj  +  A“(X  —  X*)  +  \  A~(X  -X*)2|  WjWj 
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Since  V^Wj)  =  0  at  the  critical  load,  we  expect  that  V2(Wj)  >  0  at  loads  slightly  less  than  the 
critical  value  and  that  N^fWj)  <  0  at  loads  slightly  above  the  critical  value.  Therefore,  we  conclude 
that 


A  g  <0 


49 


According  to  (49)  the  numerator  of  (47b)  is  always  negative,  but  the  denominator  of  (47b)  is  a 
homogeneous  cubic  in  Wj  and  the  sign  is  reversed  by  a  reversal  of  the  buckling  mode.  In  this  case,  an 
adjacent  state  of  equilibrium  exists  at  loads  above  (,\  >  X*)  or  below  (X  <  X*)  the  critical  value.  In 
view  of  (49)  and  (48b),  an  equilibrium  state  above  the  critical  load  is  stable  and  a  state  below  is 
unstable, 


Now,  suppose  that 


X,  =  A3  =  0 


50 


Then,  in  view  of  (37),  (40),  (41)  and  (50),  the  potential  of  (34)  and  (36)  takes  the  form: 


v  -  ( ^Aijvivj +  j  Aij«  ViWjW,  fi  AijK  NWjW.w,  •» ...)  ♦  <\ K- 


+  e  —  AjjWjWj  +...)+  O(o^) 


51 


The  underlined  term  of  (51)  dominates  if  n  f  0  and  if  c  is  sufficiently  small.  The  term  is  odd  in 
Wj  and,  therefore,  always  provides  a  negative  potential  change  at  any  load  X  f  X*.  A  condition  for 
the  existence  of  stable  states  at  noncritical  values  of  load  follows: 


AjWj  =  0 


52 


However,  the  bucking  mode  W(  is  independent  of  the  coefficients  Aj.  Therefore,  equation  (52) 
implies  generally  that 


Aj  =  0 


53 


Now,  we  accept  the  remaining  terms  indicated  in  (51)  as  our  approximation.  Also,  in  view  of 


(50), 


o  x  • 

t-  -  n  =  X  —  X* 
2 


54 
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Our  approximation  of  (51)  follows: 
,2 


v4f4(TAiivivi  +  4  Aii«ViWiW«  *  j,Aii«lWiWiW«Wl)  +  <s2  I  Aiw  o  -»*>55 


1  .-.-■ 


i 


Again,  we  require  a  stationary  potential  for  variations  of  the  displacement  V,j.  The  equations  of  J 
equilibrium  follow:  i 


e2«AijVj--c2Aii,WjWK 


Let  k  V :  denote  the  solution  of  (56).  Then,  it  follows  that 


e2K  2A;:VjV;  =  ~  C2K  A|;„W;W  V; 


56 


IJ  I  J 


IJK  J  K  I 


57 


If  the  solution  «  Vj  and  (57)  are  used  in  (55),  then  our  approximation  of  the  potential  takes  the  | 
form: 


V  =c4A4+e2(X--  X*)A2 


58, 


where  A2  is  defined  by  (44)  and 


A4  =4,  AijK|WiWjWKW|---AijViVj 


59 


The  solution  of  (56)  determines  the  unit  vector  Vj  which  renders  V  stationary,  but  still 
dependent  upon  the  distance  c.  The  principle  of  stationary  potential  gives  the  equilibrium 


condition: 

dV„ 

O 

—  x*ja2  ~  o 

J 

i 

‘ 

de 

4t-  A4  +  2e(X 

■ 

’  60 

or 

c2  ...  a2 

1 

(x  X*) 

• 

61 

2A4 

i 

According  to  (33a, b),  (49)  and  (61),  a  stable  adjacent  state  of  equilibrium  can  exist  only  at  loads 
above  the  csitical  value  (X  ^-'X*)  and  a  state  below  the  critical  value  is  unstable. 

i  1  i 
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!  Koitor  provides  rigorous  arguments  for  the  conditions  (53)  and  (49)  if  the  critical 
configuration  is  a  stable  equilibrium  configuration  for  loads  less  than  the  critical  value.  For 

1  i  I 

example,  the  two-dimensional  system  has  equilibrium  configurations  which  trace  a  line  parallel  to 
the  X  axis,  as  shown  in  Fig.  6.  The  portion  OP  represents?  stable  states,  the  bifurcation  point  Pthe 
critical  state,  PR  represents  unstable  states  ofi  the  reference  configuration  and  PQ  the  stable 
postbuckled  equilibrium  states.  Here,  the  principle  of  stationary  energy  in  the  critical  configuration 
at  any  load,  Heads  to  the  equation  (53)  and  the  principle  of  minimum  energy  in  the  stable  states  of 
OP  (Jt<  X*)  leads  to  the  inequality  (49). : 

I 

I  I 

!  >  ! 

If  the  cubic  term  of  V  does  not  vanish,  then  equilibrium  states  trace  paths  with  slope  X'at  the 

i  s  - 

critical  load,  as  shown  in  Fig.  6a.  If  the  cubic  term  vanishes,  then  X  =  0  and  the  equilibrium  states 
trace  paths  as  shown  in  Fig.  6b.  In  each  figure,  the  solid  lines  are  stable  branches  and  the  dotted 

.  *  5  5  i 

lines  are  unstable.  .  , 


>  P 


(a) 


■  i 


Fig.  6 


A3=*  0 

Aj  #  0 
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Practically  speaking,  many  structural  systems  display  the  instability  patterns  of  Fig.  6,  that  is, 
the  prebuckled  configuration  of  the  ideal  structure  is  an  equilibrium  state  under  all  loads.  Notable 
examples  are  the  column  under  axial  thrust,  the  spherical  or  cylindrical  shell  under  external  pressure 
and  the  cylinder  under  uniform  axial  compression.  Essentially,  each  retains  its  form  until  the  load 
reaches  the  critical  value  and,  then  buckles.  In  the  case  of  thin  shells,  initial  imperfections  cause 
pronounced  departures  from  the  initial  form  and  often  cause  premature  buckling  (X  v  A*). 

Our  analysis  of  stability  at  the  critical  load  is  limited.  The  reader  should  note,  especially,  that 
any  of  the  various  terms  of  the  potential,  e.g.  V2,  V4,  may  vanish  identically.  Then,  further 
investigation,  involving  terms  of  higher  degree,  is  needed. 

EFFECT  OF  IMPERFECTIONS  UPON  THE  BUCKLING  LOAD 

In  the  monumental  work  of  Koiter^,  an  important  practical  achievement  was  his  assessment  of 
the  effect  of  geometrical  imperfections  upon  the  buckling  load  of  an  actual  structure.  Here,  we 
outline  the  procedure  and  cite  the  principal  results: 


Under  the  conditions  of  dead  loading  upon  a  Hookean  structure  ,  the  energy  potential  V  of 
the  actual  structure  is  expressed  in  terms  of  a  displacement  u,  from  the  critical  state  of  the  ideal 
structure  and  a  parameter  e  which  measures  the  magnitude  of  the  initial  displacements  of  the  actual 
unloaded  structure: 


<v 

1,1  1 

1 

V  = 

2  AjjMuiuj  +  21  AjjK(X)ujUjuK  +  ^  AjjK|(A)UjUjUKU| 

+  e 

Bj(A)Uj  +--  Bjj(A)UjUj  +  ... 

62 


Here,  the  linear  terms  in  Uj  vanish  in  the  first  bracket,  because  the  reference  configuration  is  an 
equilibrium  configuration  of  the  ideal  structure  (e  =  0)  at  any  load. 


As  before,  the  components  u(  and  (X  -  X*)  are  p*  yanded  in  powers  of  the  arc  length  t  along 
the  ideal  curve  of  Fig.  7.  Here,  we  make  an  assumption  that  the  initial  deflection  of  the  actual 
structure  is  nearly  the  buckling  mode  W,  of  the  ideal  structure.  Therefore,  we  have  the 
approximation: 


.  —  ?  K 

U:  =  e W;  +  ~  V; 

2 
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o  K 

(X  -A*)  =  cA'  +  (.•  ~ii 


64 


20 


In  the  present  case,  the  second-order  terms  of  (63;  and  (64)  contain  an  unspecified  parameter  K, 

A 

because  these  terms  do  not  represent  deviations  from  the  tanger.i  (eW)  along  the  ideal  path  of  Fig. 
7,  but  represent  the  displacement  (d)  which  carries  the  system  to  the  actual  path  as  depicted  in  Fig. 
7. 


Upon  substituting  (63)  and  (64)  into  (62)  and  acknowledging  (16),  (18)  and  (20\  we  obtain 

i  Aij«  I WjWjW^W-  +  ~  AjjVjVj  +  ~  A^W.WjV, 

K 


a3 +  A^X' 

+  e4 

. 

L 

^ _ J  _ _ 

+  -  A-jW-Wtp  +-  A^W-WdX') 
4  >J  I  j  ?  'i  1  J 


+  0(e5)  +  e[,3iWi  •»  e2Bj  “Vi 


65 
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where  A3  and  A2  are  defined  by  (43)  and  (44),  as  before.  Since  the  relative  magnitudes  of  e  and  e 

O 

are  unsp°cified,  we  must  suppose  that  the  terms  0(e  )  and  0(ee)  dominate  (65),  if  A3  f  0  and  X  f  0. 
Then,  we  have  the  approximation: 

V  =  e3A3  +  e2A2(X-X*)+eeBjWj  66 

The  stationary  condition  of  equilibrium  follows: 

~  =  3e2A3  +  2eA'2(X  -  X*)  +  eBjWj  =  0  67 

Now,  recall  that  V  is  not  the  potential  increment  from  the  critical  configuration  of  the  actual 
structure  but  the  potential  referred  arbitrarily  to  the  critical  configuration  of  the  ideal  structure.  An 
equilibrium  configuration  of  the  actual  structure  is  stable  or  unstable,  respectively,  if  the  potential 
is  a  minimum  or  a  maximum;  therefore,  the  critical  load  X  of  the  actual  structure  satisfies  the 
conditions 


68 


>  0  =»  stability 

69a 

<  0  =<>  instability 

69b 

Observe  that  the  distinction  between  stability  and  instability  of  a  post-buckled  state  rests  upon  the 
same  conditions,  (69a)  and  (69b),  as  the  ideal  structure  [see  equation  (48a)]  and  that  the 
conditions  are  independent  of  the  imperfection  parameter  e. 


If  A3  f  0,  then  the  imperfect  structure  deflects  and  reaches  a  critical  state  of  equilibrium  when 
(67)  is  satisfied  and  (68)  vanishes.  If  (X  —  X*)  is  eliminated  from  the  two  equations,  then 


3e2Ao 

-m 


70 


The  sign  of  'he  sum  BjWj  is  arbitrary,  since  a  change  of  sign  is  effected  by  redefining  the  parameter 
e.  Therefore,  we  can  choose  e  so  that  B,W,  has  the  opposite  sign  of  A3.  Then  the  condition  (70)  for 
a  critical  state  is  fulfilled  only  if  e  ^  0.  From  our  observations,  we  know  that  an  imperfect  structure 
tends  to  buckle  in  a  preferred  direction,  depending  upon  the  character  of  Ihe  geometrical 
deviations.  In  the  present  case,  if  A3  <  0,  a  critical  state  occurs  only  if  e  <  0.  A  plot  of  load  versus 
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deflection  is  depicted  in  Fig.  8a;  here,  a  negative  value  e  produces  buckling  in  a  negative  mode  (cW 

A/ 

<  0)  according  to  the  curve  OP,  whereas  the  positive  value  e  produces  only  stable  states  along  the 
path  OP. 

A  real  structure  which  behaves  in  the  manner  of  Fig.  8a  is  the  frame  of  Fig.  8b.  If  the  vertical 
strut  is  bent  to  the  right  or  left  the  imperfection  parameter  e  is  negative  or  positive,  respectively. 
The  rotation  0  of  the  joint  serves  as  a  generalized  coordinate^  =  q)  and  plots  of  load  versus  rotation 
take  the  forms  of  Fig.  8a.  The  frame  under  eccentric  loading  has  been  studied  experimentally  by 
Roorda^  and  theoretically  by  Koiter®.  The  latter  computations  show  remarkable  agreement  with 
the  former  experimental  results. 

Now,  let  us  turn  to  the  structure  of  Fig.  6b,  characterized  by  the  conditions 

A3^  X'-O  71 


Fig.  8 
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Now,  terms  0(e3)  are  absent  from  the  potential  of  (65).  The  latter  must  be  stationary  with  respect 
to  the  displacement  Vj;  for  equilibrium, 

£2KA.jVj  ^WjW.-ZeBj  72 

If  KVj  denote  the  solution  of  (72),  then 

e^WjW,,  V,  "-e^AjjVjVj  -  2eKBiVi  73 


In  accordance  (64),  (71)  and  (73),  our  approximation  of  the  stationary  value  of  (65)  follows: 


V  = 


..3 


[eA^  +  A2(X  ~  X*) 


+  ecBjWj 
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As  before,  the  potential  V  is  still  dependent  upon  the  distance  e.  The  stationary  condition  of 
equilibrium  follows: 


dV 

de 


=  4  c3A4  +  3e2A2(X  -  X*)  +  eBjWj  =  0 


The  stability  of  equilibrium  depends  upon  the  second  derivative  as  follows: 


c i^V 

dc 2 


=  12c2A4  +  6t' A^(X 


X*) 


0  stability 
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76 

77a 


<  0  *=»••  instability 
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Again,  the  critical  state  of  equilibrium  is  characterized  by  vanishing  of  the  first  derivative  (75)  and 
second  derivative  (76).  The  elimination  of  (X  —  X*)  yields  the  result: 


2c3A4 

BjWj 


78 


Again,  we  note  that  the  definition  of  e  and  the  sign  of  Wj  are  arbitrary  and,  therefore,  we  assume 
that  Wj  renders  the  sum  BjWj  >  0.  Then,  the  condition  (78)  for  a  critical  load  is  attained  if  A4  •  0, 
e  <  0,  in  keeping  with  (33b).  Now,  the  structure  also  exhibits  instability  at  the  critical  load  if  the 
sign  of  the  parameter  e  and  the  buckling  mode  are  both  reversed.  A  plot  of  load  versus  deflection  is 
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depicted  in  Fig.  9.  The  linkage  of  Fig.  1  and  the  curves  of  Fig.  2  exemplify  such  structural  systems. 
The  structure  is  symmetrical  and  can  buckle  in  either  direction  depending  upon  the  character  of  the 
geometrical  deviation.  In  either  case,  it  is  stable  or  unstable  depending  upon  the  sign  of  the  constant 
A^.  The  critical  load  X  of  the  actual  structure  may  be  much  less  than  the  critical  load  X*  of  the  ideal 
structure. 

If  the  distance  e  is  eliminated  from  (67)  and  (68)  =  0,  or  if  e  is  eliminated  from  (75)  and  (76)  = 
0,  the  critical  load  X  =  X  is  expressed  in  terms  of  the  imperfection  parameter  (-e) : 

J_ 

X  =  X*  -  n(-e)(BjWj)(-An)n‘2  n4  (-A^"*1)  79 

Here  n  =  3  if  A?(  <  0,  n  =  4  if  A3  =  0  and  A4  <  0,  and  e  <  0.  A  plot  of  the  actual  buckling  load 
versus  the  imperfection  parameter  has  the  appearance  of  Fig.  10.  Since  the  curve  is  tangent  to  the 
axis  X  at  e  =  0,  small  imperfections  can  cause  considerable  reduction  of  the  buckling  load. 


Fig.  9 


THEORY  FOR  CONTINUOUS  BODIES 

1  O 

The  stationary  criteria  '  were  originally  given  for  continuous  bodies.  The  concepts  are  the 
same,  but  the  mathematical  form  is  altered:  Our  discrete  displacements  u:  are  replaced  by  the 
continuous  field  Uj(xj)  and  the  sums  of  (11)  by  integrals  of  the  continuous  field.  The  stationary 
conditions  of  Trefftz^  and  Koiter^  apply  to  the  integrals  of  corresponding  degree,  e.g.  Trefftz 
criteria  is  the  stationary  condition  upon  the  homogeneous  functional  of  second  degree  which 
replaces  the  quadratic  function  V2  of  (14). 

Koiter's  theory  has  contributed  most  significantly  to  our  understanding  of  instability, 
snap-buckling  and  imperfection  sensitivity  of  thin  shells.  The  reader  is  referred  to  the  work  of 
Koiter^;  Budiansky  and  Hutchinson®,  Hutchinson^,  Hutchinson  and  Amazigo*®  and  B. 
Budiansky^. 

Finally,  we  note  that  the  continuous  structure  can  always  be  approximated  by  a  disciote 
system  whereupon  the  foregoing  equations  are  applicable.  In  particular,  a  shell  can  bo  subdivided 
into  finite  elements,  the  deformation  can  be  approximated  by  interpolation,  and  then,  the 
stationary  criteria  can  be  applied  to  the  discrete  model  in  the  manner  of  Rayliegh-Ritz.  The  success 
of  such  methods  will  depend  greatly  on  the  complexity  of  the  continuous  body,  the  buckling  mode 
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